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^ ■ Abstract 



Using the dynamical triangulation approach we perform a numer- 



O 

Q I ical study of a supersymmetric random surface model that corre- 

O I sponds to the large N limit of the four-dimensional version of the 

■ IKKT matrix model. We show that the addition of fermionic de- 

-4— > . 

■ grees of freedom suppresses the spiky world-sheet configurations 
that are responsible for the pathological behaviour of the purely 

D I bosonic model. We observe that the distribution of the gyration 

I radius has a power-like tail p{R) ~ We check numerically 

that when the number of fermionic degrees of freedom is not susy- 
^ ■ balanced, p{R) grows with R and the model is not well-defined. 

Numerical sampling of the configurations in the tail of the dis- 
tribution shows that the bosonic degrees of freedom collapse to a 
one-dimensional tube with small transverse fluctuations. Assum- 
ing that the vertex positions can fluctuate independently within 
the tube, we give a theoretical argument which essentially ex- 
plains the behaviour of p{R) in the different cases, in particular 
predicting p{R) ~ R~^ in the supersymmetric case. Extending 
the argument to six and ten dimensions, we predict p{R) ~ -R~^ 
and p{R) ~ R~^^, respectively. 

Introduction 

During the last decades we have witnessed a rapid development of methods of 
geometrical quantization. It was triggered by the discovery of the integration 
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measure over two-dimensional geometries which was apphed in calculations 
of quantum amplitudes in string theory In practice, it turned out that 
the geometrical quantization of strings works only for either non-physical 
dimensions (i < 1 0, or in the critical dimension d = 10 where the world- 
sheet degrees of freedom of the string decouple from the theory 0, |^ . 

The case of non-critical strings with d < 1 was also independently for- 
mulated in terms of a matrix model which is equivalent to the dynamical 
triangulations approach This method allows for calculating different 
topological contributions in the double scaling limit However, both the 
Polyakov conformal field approach and the matrix model calculations break 
down at d = 1. The problems that turn up at this point, also known as the 
= 1 barrier', are attributed to instabilities of the conformal mode. 

Whether it is possible to stabilize the conformal mode for surfaces em- 
bedded in physical dimensions is not yet clear. Some attempts to solve this 
problem have been made within the conformal field formalism 0. Another 
idea is to extend dynamical triangulations to the sypersymmetric case, and 
to attack the problem directly in physical dimensions. By introducing super- 
symmetry one hopes to avoid the instabilities of bosonic surfaces that mani- 
fest as the degeneration of the world-sheet geometry into branched polymers 

i- 

Disrectization of supersymmetric surfaces is known to be a difficult prob- 
lem. A discretization of the world-sheet supersymmetry, for example, inherits 
all the difficulties which appear arleady on regular lattices and on top of 
this it introduces new problems related to the fact that the symmetry must 
be local 0. Supersymmetry is broken explicitly by the lattice, and it is not 
clear how to ascertain whether it gets restored in the continuum limit. One 
encounters similar problems for the Green-Schwarz type of surfaces where 
one has to preserve the K-symmetry, which is a local world-sheet symme- 
try [|12]. One also attempted to introduce supersymmetry into the surfaces 
models indirectly by considering models with effective actions with extrinsic 
curvature terms obtained by integrating out fermions from the supersymet- 
ric theory |TD[. Such discretized surfaces have been extensively studied (see 
e.g. ||ll|]), but it is not clear what is nature of the critical point observed in 
the discretized theory and whether it can be directly related to the original 
supersymetric theory of continuum surfaces. 

In this paper we study a model of supersymmetric surfaces that consti- 



tutes a sort of classical limit of the IKKT matrix model [0 obtained by 



sending the matrix size N to infinity ||13| , \T4^ . The attractive feature of the 



surface model is that the supersymmetry is introduced in such a way that in 
principle the model can be straightforwardly discretized. By re-writing the 
IIB string action in the Schild gauge |]r3|, one avoids the problems related 
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to the proper treatment of the local ^-symmetry in the discretization. The 
action that results from this procedure has no redundant degrees of freedom 
left, and no local symmetry. 

We investigate this model by means of the dynamical triangulation ap- 
proach. Some preliminary studies in this direction were already done in 
[T5| . We perform Monte-Carlo simulations of surfaces embedded in a 4- 



dimensional target space. We find a theoretical picture that explains the 
numerical results for d = 4, and allows us to make some predictions for the 
behaviour of the ten-dimensional version of the model. The results presented 
here give an insight into typical geometrical features of surfaces that play an 
important role in the ensemble generated by the IKKT matrix model ||13|| . 

The paper is organized as follows. First we shortly recall the model and 
describe the discretization scheme. Then we discuss singularities of the purely 
bosonic model and show that they are removed by adding fermions. For 
the model with fermions we estimate the large L behaviour of the partition 
function by assuming that it is dominated by tube-like configurations. We 
present the results of the MC simulations to show that tubes do indeed 
dominate in the ensemble. We conclude the paper with a short discussion, 
concerning in particular the relation of the results to the matrix model p!3| , 
|T6l |T7i. 



The model 



The action for JIB strings can be cast into the following form : 



S{g,X,^,^) = Jd'^^ (^^{X^,X''Y-'-^r,{X^,^} + X^ (1) 

Here Qab is the metric tensor on the world-sheet; X^^ are bosonic co-ordinates 
in the d-dimensional target space; ^ and are fc-dimensional spinors; and 
are k x k Dirac matrices. For d = 10, the dimensionality of the spinor 
representation is k = 32, and ^, are Majorana-Weyl spinors. For d = 4, 
we have k = 4, and ^, \l/ are Weyl spinors. The Poisson brackets in (|l|) are 
defined as 

{X,Y} = ^e'''daXd,Y. (2) 

Based on this action the authors of |jl3|] introduced a matrix model which 
was then interpreted as a constructive definition of superstring theory. The 
matrix model is obtained from (0) by substituting fields X'^(^i,^2) 
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and "^{C,!, ^2) —>■ ^ab hjNxN traceless Hermitean matrices, and the Poisson 
brackets by commutators {, } — "^i, ]• The resulting theory 

S = -Tr + hr^lAi^, $] - Al^ (3) 

is quantized by the Feynman integral 

Z= / dA d$ d$ e""^ (4) 



with the standard integration measure for traceless Hermitean matrices. This 
procedure defines a kind of third quantization of string theory, which is sup- 
posed to reproduce in a double scaling limit the sum over all topological and 
geometrical excitations of the world-sheet of strings [0, The original 



model 131 was defined in ten dimensions. 



Here we will discuss the semiclassical limit of 00, in which the 



model reduces to the string theory on a fixed topology with the partition 
function 

5(g,X,*,*) 



where the action is given by (|1]). We will concentrate for now on the case 
d = 4, where we can most easily perform numerical simulations. 

We use the dynamical triangulation approach to regularize the integration 
measure The partition function is then given by 

Z = Y,e'^''^N (6) 

N 



where 



Y^^J \{d'X,\{d'^>d'^,e-'-^^-^-''^^ . (7) 



Zn — 

The sum runs over equilateral, oriented triangulations T with triangles and 
n = N/2 + 2 vertices. Ct is the symmetry factor of a given T, and St is the 
discretized action on this particular triangulation. Primes on the products 
indicate that the zero modes were removed (see the following discussion). 
The fields Xj and are now located on the vertices of T. The action has 
one parameter, namely the cosmological constant A. The continuum limit is 
taken by adjusting A in the grand-canonical partition function Z{X) to its 
critical value and sending the lattice spacing a to zero. Alternatively, if one 
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Figure 1: Two neighbouring, oriented triangles (ijk) and (jin), connected 
by the hnk (ij). The 'outer' vertices k and n can be accessed by the operation 
u as defined in the text, with u^ij) = k and uj^ji) = n. 



prefers to work with the ensemble of canonical partition functions zat, one 
can take the continuum limit by simultaneously sending N to infinity and 
a to zero while holding the physical area of the worldsheet A ~ No? fixed. 
This latter approach is the one we will choose here. 

The discretized action consists of two terms, a purely bosonic part Sb 
and a fermionic part Sp- Following the suggestions in |T^, these are : 



(ijk) 

+2 {X^Xl + X%Xl + XlXl) } (8) 

{ij) 



In contrast to the action presented in here \l/ are Weyl fermions. 
The first sum runs over all triangles of the triangulation, the second over 
all links. The constant in (H) is set to B = {8\/3a^)~^. The operation cu, 
when acting on an oriented link (ij), gives that neighbour of i that comes 
after j when going counterclockwise around i. In other words, if we denote 
two neighbouring oriented triangles by (ijk) and (jin), then uj{ij) = k and 
Lj{ji) = n (see figure |I]). 

Two things should be noted about this action. Firstly, the bispinors 
ff, in the target space have a handedness, i.e. ^ = |(1 + r5)\l/ and 
^ = ^^{1 — Fs). Using the chiral representation of the Dirac matrices in the 



5 



Euclidean sector, 



' ro=f ° > r,= f V ] (10) 



a y ' V J ' ^ \0 -(To 

where a are the Pauh matrices and (Tq is the 2x2 unity matrix, we can 
parametrize the chiral bispinors ^ and \l/ in terms of two-component spinors 

^ = (0,^), ^=(t) ■ (11) 



^ ^ 

Thus, we can re-write the fermionic part of the action as 

SpiX, ^,^) = Y, - /S^o] (12) 

where f^^ = (^X^^^.^^ — Xj^^^.^ are real numbers that are antisymmetric in 

i and j, = -/j^^. 

Secondly, the action has a zero mode related to the translational invari- 
ance in the bosonic sector, —>■ X^ +5. As a remnant of the original 
supersymmetry, there is also a similar symmetry in the fermionic sector : 
V^j — ^ V'i + e, V'i V^j + e [iH- This might be not evident at first sight, but 
the change in the action is indeed zero, 

5Sf = ^ e (ifija - /° (To) + ^» [^fij^ - ftj^^o^ e + 

+ ^e(^^.a-/Sao)e = (13) 

{ij) 

because for each vertex i the sum of the /y over all its neighbours j is zero : 

//;• = 0. 

These zero modes have to be divided out in the measure. In practice, this 
can be done by inserting a delta function 5'^{Xn) and an additional product 
^n^n ^ which acts like a delta function for the Grassmann variables ipn 
- for one vertex. 

Taking all this into account, the canonical partition function for the four 
dimensional model now reads 



^ „ n— 1 n 

H C^j \{d'X,\{d'',pd''^i-i,n^n-e~'-^''^e-^^^^)^^^^^^^^ (14) 



, n-l 

zn = 
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where 



aij = ifija - fijao . (15) 

Integration over the fermions results in the appearance of a factor D{X) = 
det a[j, where a'^j is a 2 (n — 1) x 2(n — 1) matrix that is obtained from aij by 
crossing out the two rows and columns that correspond to the vertex n. The 
determinant D{X) is a non-negative function of bosonic coordinates X, a 
feature that is essential for the MC simulations. It follows from the structure 
of the matrix a'-, which has pairs of complex conjugated eigenvalues A, A*. 
(More specifically, if rji is an eigenvector to the eigenvalue A, then (T2?7* is an 
eigenvector to the eigenvalue A*.) 

The final form of the partition function is now 

„n-l 

.^(7) = E ^ / n d'X^ DiXr e-^-(^) (16) 

where for convenience we introduced an additional parameter 7. For 7 = 
we have the purely bosonic model; 7 = 2 is the case considered in |T3; and 



7 = 1 corresponds to our partition function ([T^). As we will see, the model 
with 7 = 1, for which the number of bosonic and fermionic degrees of freedom 
is susy-balanced, is the only interesting one, because only in this case is the 
partition function ([iBD well-defined. 



Singularities of the bosonic partition function 

The behaviour of the partition function for the purely bosonic model zn{0) 
has already been discussed in the context of the quantization of the Nambu- 
Goto string [^, where it was shown that this partition function is in general 
not integrable. Namely, one can explicitly construct triangulations for which 
the integral over bosonic fields is divergent. 

The source of this divergence can be intuitively explained as follows. The 
bosonic part of the action is the sum of all triangle areas squared, with the 
areas measured in the target space. Thus, as long as the areas of all triangles 
stay constant, we can change their shape without changing the action. As 
an example, let us pick an arbitrary vertex m on the triangulation. It has q 
neighbouring vertices that are linked to each other, forming a loop around 
m (see figure H (left)). Each link in this loop is the base of one of m's 
neighbouring triangles. Now let us shrink these bases to a small length r. 
(Obviously, this step will change the action, but this is a finite and local 
change that is bounded independently of r.) Next, we can move m outward 
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Figure 2: Origin of spikes : (Left) When the white points approach each 
other to a distance ~ r, then the black point can freely move in the embedding 
space inside a sphere of radius ~ 1/r. (Right) Example of the non-integrable 
configuration : the entropy factor from four spikes (black points) overcomes 
the damping effect of the entropy that comes from the white points which 
are forced to move close to each other. 

to a distance of about ~ 1/r = L, thus creating a long 'spike'. In this way, 
we can have an arbitrarily long spike that does not visibly affect the action 
if we just simultaneously make the bases small enough. 

The phase space volume for such a configuration is a factor L'^ coming 
from the point m times a factor r^''"^^'^ = L~^'^~^^'^ coming from constraining 
the q neighbours of m to a sphere of radius r around one of them. Thus, a 
spike of length L is suppressed, but only power-like, i. e. the spikes are not cut 
off at any particular scale and arbitrarily long fluctuations can occur. When 
instead of a single-vertex spike one considers configurations with multiple 
spikes, such as the one in figure (right), one can check that the spike length 
distribution may have a logarithmic or stronger divergence. 

Singularities of the susy partition function 

The question that now arises is, what can be expected to change in this 
picture if we add the fermionic part to the action and consider the complete 
partition function zn{^)'^ Intuitively, it seems that fermions should introduce 
a repulsive core to the effective potential between neighbouring vertices on the 
triangulation, and that this repulsion should prevent vertices from occupying 
the same position in the target space, which in turn should keep the loop 
around each vertex from shrinking beyond a certain minimal length. In this 
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case, to create a spike would cost energy, because now the area of a triangle 
would have to grow to allow it to become elongated. Therefore, one can 
hope that the introduction of fermions might suppress spikes, along with the 
divergences they cause. 

Introducing the fermions is equivalent to replacing the bosonic action Sb 
with an effective action 

SB,eff{X) = Sb{X) - \ogD{X) (17) 

where D{X) = det a is the determinant of the matrix defined in (|T5|). 
D{X) vanishes if for any vertex all its neighbouring vertices meet in a single 
point, because then the two rows of the matrix that correspond to this vertex 
become zero. Similarly, if all neighbouring vertices do not exactly meet in 
one point but move very close to each other, then all entries in these two rows 
become very small, and the determinant as a whole should likewise become 
small. In other words, we expect the determinant to discourage small loops 
and hence also spikes. It is natural to expect that short loop lengths will be 
suppressed by a power-like fall-off. 

To test this, we used standard Monte Carlo methods to simulate the 
model with fermions. The triangulations, weighted by the effective action 
(|T7|) , were produced by a Metropolis algorithm for geometrical flips and up- 
dates of the fields X^. In four dimensions the effective action is well-defined 
since the determinant is positive definite as discussed before. The bulk of the 
computer time is spent on calculating the determinant D{X). This allowed 
us to produce reasonably large statistics only for fairly small systems (up to 
16 vertices). 

The simplest quantity to measure is the average action. Using the fact 
that the bosonic action and the fermionic determinant are uniform polyno- 
mials of the bosonic co-ordinates, one can show by a simple scaling argument 
that the average bosonic action for (i = 4 is {Sb) = (^ — 1)(1 + 7/2). We 
used this formula as a rough test of the program. 

We also measured the gyration radius 



-V(X,-Xcm)2 (18) 
n ^-^ 

i 

and the length of the smallest loop around a vertex that appears on each 
configuration. ('Length', here, means the sum of the squared lengths of all 
links in the loop.) Figure ^ (left) shows that in the purely bosonic model 
there is a strong correlation between the appearance of short loops on the tri- 
angulation and a large gyration radius, which supports the picture discussed 
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Figure 3: The gyration radius versus the shortest loop around a vertex on 
the triangulation, for the model without fermions (left) and with fermions 
(right). The number of vertices is n = 16. 



in the previous section. Note that all results for the bosonic model should 
be taken with a generous helping of salt, since as noted above the partition 
function is actually not defined. These results are presented here merely to 
illustrate the fact that the source of the singularities can really be found in 
the appearance of local spikes (see figure ^ (left)). For the fermionic case, 
small loops are clearly suppressed, and the only correlation between loop 
lengths and gyration radius is a trivial one, namely that if the shortest link 
on a triangulation has a length r, the gyration radius itself cannot be smaller 
than r, either. We conclude that fermions do in fact suppress spikes. For the 
smallest possible spherical triangulation, the tetrahedron, one can actually 
show analytically that the partition function for the fermionic model exists 



in four dimensions, whereas for the bosonic model it does not |15]. 

However, both the analytical calculations for the tetrahedron and numer- 
ical simulations of larger surfaces indicate that even in the fermionic case the 
gyration radius distribution has a power-like tail (see the appendix and fig- 
ure ^ (left)). What is the source of the singular behaviour in this case, when 
we know that we do not have local spikes? It turns out that there is another 
type of configuration that can contribute to the tail of the probability distri- 
bution p{R) for large R. These things look like needles; the four-dimensional 
degrees of freedom of the co-ordinates in the target space collapse to a one- 
dimensional, elongated narrow tube. To indentify these configurations we 
used the principal component analysis (PCA). For each analysed configura- 
tion we calculated the correlation matrix Q^i, = ^ X'^X'^ and its eigen- 
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Figure 4: Snapshots of a susy triangulation from the tail of the gyration 
radius distribution (right), and a bosonic configuration of comparable size 
(left). The x-component in the plot corresponds to the projection of the 
vertex position onto the longest principal axis in the system, while the y- 
component describes the distance from the axis. L denotes the length of the 
tube, r its radius and R the gjnration radius 



values. The eigenvector corresponding to the largest eigenvalue defined the 
principial axis. For each point on the configuration Xi we then calculated its 
projection on the pricipial axis, X^, and its distance from it, di. We defined 
the length L and radius r of the tube by 

L = maxXj^ — minX^ and r = max di (19) 

i i i 

The ratio of length to width can serve as an estimator of the type of the 
configuration. For the bosonic string, this ratio was 2.26(8), and for the 
fermionic one, 212(20). Two typical results of this analysis are depicted in 
figure ^ where we show, for comparison, snapshots of triangulations for the 
bosonic (left) and susy (right) models for n = 16, projected onto the principal 
axis. 

How can the appearance of needles be explained? The bosonic part of the 
action (|]) is a sum of contributions from separate triangles. The contribution 
from a single triangle is minimal when its three vertices lie on a line, in which 
case the square of the triangle area is zero. One type of zero-action triangles 
are the spikes that were discussed in the previous section. On a spike, two 
vertices of the triangle lie very close to - or, in the limiting case, on top of 
- each other, so that together with the third vertex they will automatically 
form a line. Another possible type of zero-action triangles are those that 
have three distinct (and distant) points that lie on, or close to, a single line. 
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Prom the point of view of just one triangle, this type of zero-action configu- 
ration is more hkely to occur than a spiky one, for which one has to move two 
vertices into exactly the same point; but when one looks at the triangulation 
as a whole, one sees that the probability for this to happen is actually much 
smaller, because having three long links means that either there will be three 
neighbouring triangles with large areas (which is suppressed by the action), 
or there will be three more triangles lying on exactly the same line. The 
same argument can then be applied to the neighbours of these triangles as 
well, and on and on until one arrives at a configuration where all triangles 
lie on a single line, i.e. a configuration that is just one-dimensional. Con- 
trary to spikes, which are defects that involve only a few triangles and that 
are basically independent of the rest of the triangulation, a needle requires 
a global arrangement of vertices that minimizes the action. Such a global 
arrangement is entropically disfavoured and therefore does not occur in the 
bosonic model. Here, however, the fermionic determinant provides an addi- 
tional factor that favours configurations with vertices far removed from each 
other. As it turns out, this term is sufficient to exceed the damping effect 
coming from the entropy. 

To see this, we will try to estimate the contribution of the needle-like 
configurations to the partition function. Let us assume that we fix one point 
on the surface (thus getting rid of the zero mode), and also set the length 
of one link emerging from this point to L. In the following, we will refer to 
the direction along this link as long and to the transverse directions as short. 
Consider now a tube of radius r = 1/L around the long direction. Clearly, 
if we constrain all the points to lie inside the tube, the area of each triangle 
will not exceed \Lr — \. As the action is the sum of the areas squared of 
the individual triangles, it will be bounded from above by \nT independently 
of L. Assuming that all the vertices are free to move independently inside 
the tube, we can estimate the contribution to the partition function coming 
from this configuration as : 

Zt^^ JdLL^ (dL dr^)""' D{L, r)^. (20) 

Here, D{L,r) is the fermionic determinant for the tube-like configuration. 
We have to integrate over only n — 1 vertex positions because we already 
fixed one of them. By choosing a convenient co-ordinate system, we can set 
this fixed vertex to the origin, and the vertex at the other end of the link 
to (L, 0,0,0). This naturally defines the direction of the needle. Since the 
needle may point in any direction in four- dimensional space, integration over 
the position of this point contributes a factor dLL^ to the integral. Each 
of the n — 2 remaining vertices then contributes one integration over L and 
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three integrations over short components, resulting finally in pO|). 

The determinant -D(L, r) is a uniform polynomial of degree 2{n—l), so one 
would expect at first glance that in the large L limit the leading contribution 
should be L'^^"'~^\ However, as we will show in the appendix, because of 
additional symmetries of the matrix all terms containing only long directions 
cancel, and we get D{L,r) ~ (r^L""'^)^ for even n and D{L,r) ~ (r^L"~^)^ 
for odd n. Inserting this into (|20|) and integrating over all but one long 
component - which corresponds to substituting J dL ^ L and J dr 1/L 
- we eventually obtain : 

{a = 2n{l — 7) + IO7 — 7 for even n 
(21) 
a = 2n{l — 7) + I87 — 7 for odd n 

Because the configuration is essentially one-dimensional and the length of 
the tube is given by L, we expect the gyration radius to be proportional to 
this, R ^ L. Thus, for large i?, we expect a similarly power-like tail in the 
gyration radius distribution, 

di? p{R) ~ di? i?-" (22) 

z/the needle-like configurations dominate in the partition function. 

For the purely bosonic case 7 = 0, a = 2r;, — 7 grows with the size of the 
system and the tail becomes less singular. This suppresses needles. In fact, 
as we already discussed, in this case the spiky configurations win entropically 
over the needles, and it is they who are responsible for the singularities of 
the partition function for bosonic surfaces. 

For 7 = 1, the formula gives a = 3 for configurations with an even number 
of vertices and a = 11 for those with an odd number. These exponents do 
not grow with n. 

For 7 > 1, the exponent a decreases with n and very quickly becomes 
smaller than one. The tail of p{R) becomes non-integrable, signalling 
that the partition does not exist in this case except for a very few small 
systems that have a larger than one, such as the tetrahedron, for which the 
formula predicts a = 1 + 27. (This result for the tetrahedron can actually be 
derived in a straightforward analytic calculation, as shown in the appendix.) 

To test these predictions, we measured the exponent a numerically for 
different values of n and 7. Numerical simulations for problems with a power- 
like distribution are difficult because in this case one has rare but important 
events. In our case, the distribution p{R) ~ -R~" means that we expect 
rare, long excursions in the configuration space which produce large i?'s. 
Technically, it is difficult to explore this part of the phase space, because the 
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Figure 5: The distribution p{R) for n = 16 and 7 = 1 (left) and n = 8 and 
7 = 2 (right). 



bulk of the probability distribution is concentrated around small R where 
the program spends almost the whole time. To prevent this, we introduced 
a lower hmit Rmin on R, as well as an upper limit Rmax to prevent the 
program from doing too long excursions to large values of R. Otherwise, the 
program would not have a well defined autocorrelation time. Experimentally 
we found it convenient to set the limits Rmin = 4 and Rmax = 10 when we 
expected a negative a, and only an upper limit Rmax = 7 when we expected 
the distribution p{R) to grow with R. In figure ^ we show two examples of 
the numerically obtained distribution p{R) : for n = 16 and 7 = 1, where 
according to the formula (|21| ) we expect a = 3, and for n = 8 and 7 = 2, 
where we expect a = —3. (Note that in the latter case, the partition function 
would be ill-defined without the upper limit imposed on R. We consider it 
here only for the purpose of testing our formula.) The numerical results for 
a are summarized in table |I]. 

We can divide these results into three different cases. For 7 = 2, where 
needles are expected to not just create a tail in the distribution but instead 
dominate it entirely, the agreement is wonderful, and there is no question 
that our formula is correct. For 7 = 1 and configurations with an even num- 
ber of vertices, the results are somewhat off, but still in the ballpark of what 
we expect. One possibility is that this discrepancy is due to finite size effects, 
and would diminish if we were to go to larger systems; another is that our 
theoretical argument is simply too rough and works only as a first approxi- 
mation. In any case, we know from the analysis of sample configurations, as 
presented in figure ^ that we do have needles in this case. For 7 = 1 and an 
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Table 1: Values of the exponent a for various combinations of n and 7. 

odd number of vertices, the difference is huge but not surprising, since the 
analysis shows that these configurations are not needles, but look rather like 
small clumps with a single long spike growing out of them. It would seem 
that these things are suppressed by a power that is smaller than —3 but still 
larger than —11, so that they can dominate over needles in odd-ra cases but 
not in even-n ones. 

The estimate of the large R behaviour of the gyration radius distribution 
for needle-like configurations can easily be extended to the d = Q and ci = 10 
cases. Since we saw in the four-dimensional case that needles are dominant 
only for even n, we restrict ourselves to these systems. If we repeat the 



counting of long and short degrees of freedom as in (pO]) , we get with r ~ 1 /L : 



■'tube 



J dL U'-^ {dLdr^^-^^f^ D{L,ry ^ j dL 1'+^'-^^^''-^^ 0(1,1/ 

(23) 



where D{L,r) is now the determinant or the Pfaffian of the Weyl fermion 
matrix in d = 6 and d = 10, respectively. The function D{L,r) is a uniform 
polynomial of the bosonic components of order {d — 2){n — 1) for d = 4,Q, 10. 
Because of zero modes occurring in the minimal blocks of the fermionic ma- 
trix (as discussed in detail for d = 4m the appendix), the terms with maximal 
power of L cancel and the leading terms for large L are (r^L"~^) ^. Thus, 
for the tube we get D(L, 1/L) ~ /,{"~5)('^~2)_ Inserting this into (p3D leaves 
us with 

Zi,, ^ j dL ^l-2(d-2)+(,-l)(n-5)(d-2) _ (24) 

In the susy case, i. e. 7 = 1, the tail of the probability distribution p{L) ~ 
^i~2{(i-2) jg intiependent of n. It is the same behaviour as the one obtained 
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in the matrix model from the analysis of the spectrum of large eigenvalues 

For 7 > 1, the needle-like configurations lead to a non-integrable singu- 
larity and the partition function is divergent in this case. For 7 = 0, the 
singularity becomes softer when n goes to infinity. We know, however, that 
in this case needles do not play any important role in the picture. 

To summarize, the discussion presented above allows predictions for the 
singular behaviour of the partition function in those cases when needle-like 
configurations dominate in the ensemble. 



Discussion 

We have investigated geometrical properties of the large N limit of the 
IKKT matrix model, which corresponds also to the reduced supersymmetric 
Yang-Mills theory We showed that the originally four-dimensional the- 
ory reduces to a one-dimensional one and is dominated by elongated tube- 
like configurations. The gyration radius distribution has a power-like tail 
p{R) ~ -R~", where a is roughly consistent with the value of 3 obtained from 
the theoretical power-counting for needle-like configurations. Repeated in 
d = 6,10 the power-counting for such one-dimensional configurations gives 
distributions p{R) ~ and p{R) ~ R~^^, respectively. Strikingly, the 
same behaviour is expected from the analysis of the large R singularities 
of the matrix model Yang-Mills integrals This result has never been 

proven for the matrix model (except for N = 2), but if it holds, in view of 
our geometrical picture it would mean that the Yang-Mills integrals are also 
dominated by one-dimensional structures for large R. 

It is clear from our considerations that if too many fermionic degrees of 
freedom were added to the model (7 > 1), the power a would decrease with 
n and the partition function would soon become divergent. On the other 
hand, if there were too few degrees of freedom, as in the purely bosonic 
case, the singularity of the needle-like configurations would become softer 
for larger n similarly as in the purely bosonic matrix model pO|]. We know. 



however, that in this case the surface model has a different type of singularity, 
corresponding to the local creation of spikes, caused by a sort of local flat 
direction in the action. The effect introduces a stronger singularity that, 
again, makes the partition function ill-defined. In summary, only the proper 
susy balance between fermionic and bosonic degrees of freedom guarantees 
the existence of the surface theory. 

The results presented in this work are based on numerical simulations of 
rather small systems. It would be very interesting to perform a systematic 
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analysis of the finite size effects. This would, however, require a considerable 
effort to improve the fermionic algorithm. 
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Appendix 



In the first part of the appendix, we will shortly repeat calculations of the 
partition function 2:4(7) for the tetrahedron |jl5|, and extend them to calculate 
the probability distribution of the link length. 

Introduce a co-ordinate system in the target space such that the origin co- 
incides with one of these vertices, setting e. g. X4 = (0, 0, 0, 0). Next, choose 
the four axes - called here L,x,y, and z - such that one other vertex lies on 
the L axis, one lies in the (L, x) plane, and one in the (L, x, y) hyperplane. 
Thus, we have Xi = (Li, 0, 0, 0), X2 = (L2, X2, 0, 0), and X3 = (L3, X3, 1/3, 0). 

The determinant D{X) in this co-ordinate system becomes fTl 



D{X) = Ll xl yl 



(25) 



Inserting this into (|T6D and integrating out all co-ordinates except Li and X2 
leads to 



^4(7) 



00 noo 

dLi / dx2 
Jo 



( 



^1-^2 



\ 



\\l\Ll + xi j 



(26) 



For large Li, X2 is of order 1/Li, fixed by the combination Llx2 in the 
exponent. The integration over X2 therefore effectively corresponds to setting 
X2 ~ (including in the measure dx2)- Thus, for large Li, the resulting 
integral is dominated by the contribution 



^4(7) 



00 PCO 

dLi p{Li) ~ / dLi 
Jo 



1-27 



(27) 



The component Li of the vertex Xi is the length L of the link between 
vertices 1 and 4. Since there is nothing to distinguish between different 
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links, p{L) ~ L^^^^'' can be viewed as the probability distribution of the link 
length for large L. 

As can be seen from (PB]), for the tetrahedron the determinant contains 
short components to the fourth power, namely We will show that this 

feature of the determinant is not particular to the tetrahedron, but appears 
in any needle-like configuration with an even number of vertices, whereas for 
any needle with an odd number of vertices we find eight short components. 
The reason for this lies in the Dirac structure of the fermionic matrix. 
To see this, let us define a matrix 



A' = A',- teA[ = ( ° )-^e( ^„ ^'t ) (28) 



rj \r-ir -r 

which for e = 1 is equal to the fermionic matrix. A' = —a'. The matrices /'^ 
are antisymmetric — 1) x (n — 1) matrices, f"^ = —{f'^Y'- By construc- 
tion, for needle-like configurations all elements of f'^ are linear in the long 
components L, while those of the f'^ are linear in the short ones, r^. Our aim 
is to prove that the leading term of the determinant det A' is proportional 
to at least for even configurations and for odd, which is equivalent to 
saying that the leading term contains at least four/eight small components, 
which in turn gives D{L,r) ~ r^L^"~^ and D{L,r) ~ r^L^""^*^ for even and 
odd configurations, respectively. We will use first order perturbation theory 
to show this. 

First of all, note that /'° has at least one zero eigenvalue, because for 
strictly one- dimensional configurations with Xi = (X°, 0, 0, 0) the action (|^) 
has an additional zero mode coming from invar iance under a change ipi —>■ 
tpi + eXf. (This is a discrete remnant of a more general invariance 
'ip{^) + ^g{X^{C,)) of the continuous action, with an arbitrary function (^(X^).) 
For the truncated matrix /'°, this corresponds to a zero eigenvector 



Wi 



i = l,...,n-l (29) 



Thus, we also have detAg = (det /'°)^ = 0, which already excludes terms of 
order e". 

Consider now separately the case of even n. Because /'° has the zero 
eigenvector w, A'q has at least two eigenvectors u and d, which we can write 
as 

Under a small perturbation —ieA[ the eigenvalues A of A'^ are changed to 
(A — ieA, A -|- ieA*). Denote the first order correction to A = by Aq. The 
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determinant to this order is then proportional to e^AoAgJ^'A^, where the 
product runs over nonzero eigenvalues of f^. Aq fulfills the equation 

det f ""1^^" .T^f'' ) = (31) 

\^ (fA[u d^A[d - Ao y ^ ^ 

which gives Aq = due to the antisymmetry of the /'^ : f'^w = for any 
k. Thus, the eigenvalue A = remains intact to first order, which means that 
all terms of order likewise vanish in the determinant of A'. Therefore, the 
leading terms must be at least of order e^. One can check that, generically, 
they are indeed of this order. 

For odd n, /'° is an even by even antisymmetric matrix and, as such, 
cannot have just one zero eigenvector but must have at least two, which 
we will call v and w. This in turn means that A'^ has at least four such 
eigenvectors, which we write as 



Repeating the perturbation analysis in e for the matrix A' (^), we find for 
the first order correction Aq : 



i=l 



. (33) 



We will now argue that all three terms in the sum vanish, and therefore 
Aq = 0. This follows from the peculiar structure of the matrices f'^. As 
we know, their elements are linear in the components X^, and moreover all 
the matrices have the same structure. Therefore, we can introduce a linear 
function /' such that f'^ = f{X^). Also, the eigenvector w is likewise linear 
in X^, and we can again write w'^ = w{X^). Now we have 

v^rw = v"'{X^)f\X^)w{X^) 

= v^{X^) f\X') (w{X^) + w{X')^ 

= v^iX') (r(X°) + /'(XO) [wiX') + wiX^)) 
= t;^(X°)/'(X° + X*)w;(X° + X') 

= (34) 

(In the first step, we just added f'{X^) w(X*) = 0; in the second, we used 
the fact that v'^{Xq) /'(Xq) = 0; and in the third, we used the linearity of /' 
and w.) 
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Thus, we have shown that there are still four zero eigenvalues even to first 
order of e. As with even n, we checked numerically that they are non-zero 
to the second order. Thus, the leading term in this case behaves as e^. 
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